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Abstract 


£-V0 

Although dry/moist potential vorticity ( -) is a useful physical quantity for 

P 

meteorological analysis, it cannot be applied to the analysis of 2D simulations. A 


P xV0 

convective vorticity vector (CW) is introduced in this study to analyze 2D 

P 


cloud-resolving simulation data associated with 2D tropical convection. The cloud model 
is forced by the vertical velocity, zonal wind, horizontal advection, and sea surface 
temperature obtained from the TOGA CO ARE, and is integrated for a selected 10-day 
period. The CW has zonal and vertical components in the 2D x-z frame. Analysis of 
zonally-averaged and mass-integrated quantities shows that the correlation coefficient 
between the vertical component of the CW and the sum of die cloud hydrometeor 
mixing ratios is 0.81, whereas the correlation coefficient between the zonal component 
and the sum of the mixin g ratios is only 0.18. This indicates that the vertical component 
of the CW is closely associated with tropical convection. The tendency equation for the 
vertical component of the CW is derived and the zonally-averaged and mass-integrated 
tendency budgets are analyzed. The tendency of the vertical component of the CW is 
determined by the interaction between the vorticity and the zonal gradient of cloud 
heating. The results demonstrate that the vertical component of the CW is a cloud-linked 
parameter and can be used to study tropical convection. 
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1. Introduction 


Because it is conserved in frictioniess, adiabatic flow in a dry atmosphere, potential 
vorticity is one of the most important dynamic/thermodynamic parameters. It has been 
studied to enhance the understanding of the genesis and development of weather systems 
for more than six decades since it was first introduced by Ertel (1942). However, it is not 
conserved when clouds develop and release latent heat. Moist potential vorticity is thus 
introduced by replacing potential temperature with the equivalent potential temperature, 
which is conserved in frictioniess moist adiabatic processes. Many studies have 
contributed to understanding the roles of dry and moist potential vorticity in the genesis 
and development of weather systems (e.g., Bennetts and Hoskins 1979; E man uel 1979; 
Danielsen andHipskind 1980; Thorpe 1985; Hoskins and Beirisford 1988; Xu 1992; 
Montgomery and Farrell 1993; Cao and Cbo 1995; Cho and Cao 1998; Gao et al. 2002). 

However, this important physical parameter cannot be applied to the analysis of the 
two dimensional (2D) simulation data. Dry/moist potential vorticity can be expressed as 

|*V0 

, where § is the absolute vorticity, Q the potential temperature in dry air and 

P 

equivalent potential temperature in moist air, p the air density, and V the 3D gradient 
operator. For 2D x-z flows (e.g., Tao and Simpson 1993, Wu et al. 1998, Li et al. 1999), 

^ = ( + 2Qcos cp)j + 2Qsin<p& and W6 = ^-i + ^-k , where u and w are the 

dz dx dx dz 

zonal and vertical wind components respectively, x and z the zonal and vertical 
coordinates respectively, i and k the unit vectors in the zonal and vertical coordinates 
respectively (j = kxi). Q is the angular speed of the Earth’s rotation, and cp is the 


latitude. — — — = . For 2D, equatorial flows (cp = 0), — = 0 . This 

p p dz p 

demonstrates that the 2D flows do not contribute to the dry/moist potential vorticity. The 
vertical component of planetary vorticity is ignored henceforth. 

The new vorticity vector in the 2D x-z frame is: 


p dz dx 


( 1 ) 


du dw 

where Q =■ — --- — and 6 e is the equivalent potential temperature, 2Qcos<p is much 
dz dx 

smaller than £ as indicated in Fig. 4a and is neglected in this study. The new vector 

£ x V0 e t dd t d6 

( ) has zonal (P x = =- — -) and vertical (P z = -= — -) components in the 2D x- 

p p dz p dx 

z frame. This new vector will be used to analyze 2D tropical convection based on hourly 
data from a cloud-resolving simulation. It will be demonstrated that this is an important 
vector whose variation is closely associated with that of tropical convection. The model, 
forcing, and experiment are briefly described in the next section. In section 3, 2D 
modeling data will be used to analyze the vorticity and equivalent potential temperature 

P x V0 

gradients, calculate — , derive its tendency equation, and examine the dominant 

P 

processes responsible for its variation in the tropical, deep convective regime. The 
summary is given in section 4. 


2. Model 
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The cloud-resolving model used in this study was originally developed by Soong and 
Ogura (1980), Soong and Tao (1980), and Tao and Simpson (1993). The 2D-version of 
the model used by Sui et al. (1994, 1998) and further modified by Li et al. (1999) is what 
is used in this study. The governing equations and model setup can be found in Li et al. 
(1999, 2002c). Several 2D cloud-resolving models have successfully simulated 
atmospheric thermodynamic profiles, cloud properties, and precipitation in the tropics 
during the Global Atmospheric Research Program Atlantic Tropical Experiment (GATE) 
(e.g., Xu and Randall 1996; Grabowski et al. 1996) and Tropical Ocean Global 
Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) (e.g., 
Wu et al. 1998; Li et al. 1999, 2002a, b, c; Li 2004). 

The model is forced by zonaUy-uniform vertical velocity, zonal wind, and thermal 
and moisture advection based on 6-hourly TOGA COARE observations within the 
Intensive Flux Array (IFA) region (Zhang, personal communication, 1999). The 
calculations are based on a constrained, variational method applied to column-integrated 
budgets of mass, heat, moisture and momentum as proposed by Zhang and Lin (1997). 
Hourly sea surface temperature (SST) at the Improved Meteorological (IMET) surface 
mooring buoy (1.75°S, 156°E) (Weller and Anderson 1996) is also imposed in the model. 
The model is integrated from 0400 LST 19 December 1992 to 0400 LST 29 December 
1992 (10 days total). Figure 1 shows the time evolution of the vertical distribution of the 
large-scale, atmospheric vertical velocity and zonal wind and the time series of the SST 
during the 10-day period. In this model setup, the horizontal boundary is periodic. The 
horizontal domain is 768 km with a horizontal grid resolution of 1.5 km. The vertical grid 


resolution ranges from about 200 m near the surface to about 1 km near 100 mb. The time 
step is 12 s. 

3. Results 

Figure 2 shows the time evolution and zonal distribution of surface rain rates 
during the 10-day integration. Convection occurs after a half-day of integration and 
begins to organize after another half-day of integration. Surface rain rates exceed 25 mm 
h' 1 during the integration. Convection moves rather slowly over the first 3 days due to 
minimal mean-layer zonal flow (Fig. lb); it later propagates eastward as westerly winds 
intensify. 

Figure 3 shows zonal-vertical cross sections of streamlines and the sum of the 
mixing ratios of cloud hydrometeors at 1800 and 2100 LST 21 and 0000 LST 22 
December 1992. In the zonal domain between 200 and 768 km, there are three, strong 
convective systems at 1800 LST 21 December associated with surface rainfall. The 
convection centered around 600 km is the strongest with hydrometeor mixing ratios of 
more than 5 g kg" 1 above the melting level. The streamlines show the circulations with 
upward motion are associated with the convection. Westerly winds prevail across the 
model domain below 600 mb. The convection between 200 and 350 km has no 
significant propagation and later weakens. The convective systems between 400 and 500 
km and 600 and 650 km move very little in 6 hours. However, clouds between these two 
convective systems develop by 2100 LST 21 December. The two convective systems 
merge, and the convection centered at 600 km weakens by 0000 LST 22 December. 



The vorticity £ and the zonal and vertical gradients of the equivalent potential 
temperature 6 e at 2100 LST 21 December 1992 are plotted in Fig. 4. Large variations in 

the vorticity and zonal and vertical gradients of equivalent potential temperature are 
associated with large magnitudes of hydrometeor mixing ratios, though the magnitudes of 
the zonal gradient are larger than those of the vertical gradient (note the difference in 
scale for the zonal and vertical gradients in Figs. 4b and 4c). Thus, both zonal (P x ) and 


£xV0 e 

vertical (P) components of ( ) will be analyzed. As seen in the following 

P 

analysis, the vector is closely associated with the variation of tropical clouds and is 
referred to as the convective vorticity vector (CW). 

Zonally-averaged and mass-integrated P x , P z , and sum of the cloud hydrometeor 

mixing ratios [i.e., q c ,q r , q x , q s , q g , the mixing ratios of cloud water (small cloud 

droplets), raindrops, cloud ice (small ice crystals), snow (density 0.1 g cm' 3 ), and graupel 
(density 0.4 g cm' 3 ), respectively] were calculated to examine their relationship. They are 
denoted by [PJ, [PJ, and [q c +q r + q t + q s + q g ], respectively (where 



pF]dz ) and are shown in Fig. 5. The time evolution of [P z ] is in phase 


with that of [ q c + q r + q x + q s + q ] as indicated by their correlation coefficient of 0.8 1, 
though more fluctuations occur in [PJ than in [q c + q r + q t + q s + q g ]. [PJ, however, 
does not follow the time evolution of [ q c + q r + q t + q s + q g \ and their correlation 
coefficient is only 0.18. Thus, the tendency of the vertical component of CW is closely 
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associated with the variation of tropical convection. Its tendency will be analyzed in the 
following discussion. 

To analyze the physical processes responsible for the variation of the vertical 
component of the CW, the tendency equation for P z is derived in Appendix, which is 
expressed by 

dP z . dP 2 dP z . t ,du dd e dwdd e . 1 dO e dB 

dt dx dz p dx dx dx dz p dx dx c p dx T 


Term A 


TermB 


TermM 


j £ j (M-). 

c p p dx T 


Term R 


( 2 ) 


In (2), the tendency of P 2 is determined by five terms on the right-hand side. Term A 
includes the first two terms and is related to dynamic and thermodynamic processes 
including the zonal and vertical advection of P z and the interaction between vorticity and 
the gradients of wind and equivalent potential temperature. Term B is the third term and 
is associated with the buoyancy force. Term M is the fourth term and is associated with 
cloud microphysical processes. Term R is the fifth term and includes the solar and 
infrared radiative forcing. 

The zonal mean and mass integration of (2) yields 
3[PJ _ r £ ,du dd e , cbv 30 ev , , r 1 dd e dB ] L/£ 3 ( P 18 0 g)] 


dt 


-[«(■ 


L )M- 


p 3x dx dx dz p dx dx c p dx T 


1 £^)i. 


c„ p dx T 


(3) 


The zonal mean and mass mtegration of the first term in (2) is much smaller than the 


R 


other terms in (3) and is thus excluded in (3). 

To examine the processes responsible for the variation of [P z ], time series of the 
four terms in the tendency equation (3) are plotted in Fig. 6. [Term M] follows the 
tendency of [jP] closely in both phase and amplitude with a correlation coefficient of 

0.97 and a root-mean-square (RMS) difference of 1.2x1 O' 6 Ks' 2 , which is much smaller 
than the standard deviation of the tendency of [i^] (3.7 xl0~ 6 Ks 2 ). [Term B] is positive 
and has moderate amplitudes, whereas [Term A] is negative. Their correlation 
coefficients with the tendency of [P z ] are smaller than 0.2, and their RMS differences 
with the tendency of [P z ~\ are similar to the standard deviation of the tendency of [PJ. 
[Term R] is negligible. Thus, 

(4) 

dt c p p dx T 

The dominance of [Term M] in the variation of[F x ] indicates that the variation of 
the vertical component of the CW is controlled by the interaction between the vorticity 
and zonal gradient of cloud heating. Eq. (A4b) shows that the cloud heating comes from 
the ice microphysical processes including vapor deposition, evaporation of liquid water 
from the surface of graupel, accretion of cloud water by snow, and melting of snow and 
graupel. This implies that Pig could be zero in water clouds. Thus, the variation of the 
vertical component of the CW is associated with the variation of the mixing ratios of ice 
clouds through the interaction between the dynamics and ice microphysics. 

4. Summary 


Q 



A new vector (convective vorticity vector, CW) is introduced to study 

p 

tropical convection. In the 2D x-z frame, the CW has two components: a zonal 
component associated with the product of the vorticity and the vertical gradient of 
equivalent potential temperature and a vertical component related to the product of the 
vorticity and the zonal gradient of equivalent potential vorticity. Both components are 
calculated using hourly data from a 10-day, 2D cloud-resolving simulation, which is 
forced by the vertical velocity, zonal wind, horizontal advection and sea surface 
temperature observed and derived from the TOGA CO ARE. 

The analysis of zonally-averaged and mass-integrated quantities shows that the 
variation of the vertical component of the CW follows the variation of the sum of the 
mixing ratios of cloud hydrometeors with a correlation coefficient of 0.8 1, whereas the 
correlation coefficient between the zonal component of the CW and the sum of the 
mixing ratios is 0.18. This indicates that the vertical component is a good indicator for 
the variation of 2D tropical convection. The tendency equation for the vertical component 
of the CW is derived to determine the dominant physical processes responsible for its 
variation. Analysis of the zonally-averaged and mass-integrated tendency budget reveals 
that the interaction between the vorticity and zonal gradient of cloud heating accounts for 
the tendency of the vertical component of the CW, indicating the interaction between 
dynamics and cloud microphysics links the CW and clouds. The analysis further shows 
that the ice microphysics dominates the cloud heating, implying there are dynamic 
differences in CW between water and ice clouds. Processes associated with buoyancy 
forcing and dynamic and thermodynamic interaction have a moderate impact on the 
tendency, whereas radiative processes have no significant effect. The prehminary, 2D 


in 



results suggest that the vertical component of the CW is a cloud-related parameter that 
can be used to study the interaction among dynamic, thermodynamic, cloud 
microphysical processes during the development of 2D tropical convection. 
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Appendix 


Derivation, of tendency equation of P_ 


The zonal and vertical momentum equations can be expressed by 


du du du d(0n) 

— = —u w c n , and 

dt dx dz dx 


(Ala) 


dw dw dw d(0Jt) n 

- — = -u w c — — — - + B, 

dt dx dz dz 


(Alb) 


where it = (pi p o ) K with K - Rlc p ; R is the gas constant; c p is the specific heat of dry 

0 ' 

air at constant pressure p, and p o =T000 mb; B - g{ — - + 0.61g v ~q c -q r - q t -q s — q g ) is 

°b 

the buoyancy force, and q v is the specific humidity. An overbar ( ) denotes a zonal 

mean, and subscript b denotes an initial value, which does not vary with time. A prime ( ) 
is the perturbation from the zonal mean. The dissipation terms are excluded in (Al) for 
simplicity. 

d d 

Taking — -(Ala) -(Alb) and using mass continuity equation, the tendency 

dz dx 

equation for the vorticity (man elastic approximation as is used in the 2D cloud- 


resolving model is expressed by 

, d Z d Z 1 dB 

dt p dx p dz p p dx 


(A2) 


Here p is the mean air density and is a function of height only. 

The tendency equations for potential temperature and specific humidity are 


dd d6 dd D 

— = ~u w- 1 --= eL + 

dt dx dz 7tc„ Tic 


, and 


(A3 a) 


1? 


(A3b) 


dq v u dq v dq v 

= —u w o , 

dt dx dz 


where 


Qcn ~ Ly C^CWD ~ PrEVP ) + As i^DEP + 0 A )^SDEP (T < A ) + 0 A )^GDEP (T < A ) 

~ Pmlts(P > T a ) ~ Pmltg (T > To)} + Tf {P$acw (T < A) + P$ F w (T < P 0 ) + Pgacw^P < A> ) 

+ P/ACR {T <T b ) + P gacr {T < 

A ) + PsACR ( 7<T.) + P GFR (T < TJ-P^iT > T a ) 

~ PsMLT (T>T a ) 

Pgmlt (T > T a ) + P moM (T <TJ- P^iT > T m ) + P mw (T oa <T< TJ}, 

(A3c) 

and 

S v = PcND ~ PrEVP + PpEP + 0 ~ A )PsDEp(P < To) + 0 _ ^l)PcD£p(P < A) 

- Pmlts (T>T 0 )~ P mTC (T > T a ). (A3d) 

Here, T is air temperature, and T 0 — 35°C. L v , L s , and L f are latent heat of vaporization, 
sublimation, and fusion at 0°C, respectively, and L s = L v + L f . Q R is the radiative heating 

rate due to convergence of the net flux of solar and infrared radiative fluxes. Pcnd is the 
growth rate of cloud water by the condensation of supersaturated vapor. Pdep, Psdep, 
Pgdep are the growth rates of cloud ice, snow and graupel by the deposition of 
supersaturated vapor, respectively. Prevp is the growth rate of vapor by the evaporation of 
rain. Pmlts and Pmltg are the growth rates of vapor by the evaporation of melting snow 
and liquid water from the surface of graupel, respectively. Pjmlt is the growth rate of 
cloud water by the melting of cloud ice. Pgmlt, Psmlt, and Pbacs are the growth rates of 
rain by the melting of snow and graupel and accretion of snow, respectively. Pmw, Piacr, 
and P thom are the growth rates of cloud ice by the deposition of cloud water, accretion of 


n 


rain, and homogeneous freezing of cloud water, respectively. Psacw , Psfw, and P sack are 
the growth rates of snow by the accretion and deposition of cloud water and accretion of 
rain, respectively. Pgacr, Pgacw, and Pgfr axe the growth rates of graupel by the accretion 
of rain and cloud water and freezing of rain, respectively. The parameterized schemes for 
these cloud microphysical processes are from Rutledge and Hobbs (1983, 1984), Lin 
(1983), Tao et al. (1989), and Krueger et al. (1995), which are summarized in Tao and 
Simpson (1993) and Li et al. (1999, 2002c). 

Taking — x (A3 a) + — — x (A3b) and defining the equivalent potential 
6 c p T 


L q 

temperature 6 e as 8 exp( v v ) yields a tendency equation for equivalent potential 

Cp T 

temperature given by 


8t 


86, 86, Lf P,j6 x hr*'. 

= —u u - w — + —L >»_£. ^ c - 


8x 


dz 


1 OA 


c , T 


(A4) 


where 

^8 = Pdep + (1 - <5; i)P sdbp (T < T a ) + (1 - 6 x )P gdep (T < T 0 ) 

~ Pmlts (T> 

To ) Tmltg (T> 

T 0 ) + Psacw ( 7<t 0 )+p sfw (t<t 0 ) + Pgacw (T < TJ 
+ Piacr (T < T 0 ) + Pgacr (T < T 0 ) + Pracr (2 <T o ) + P GFR (T < T 0 ) - P^^ (T > T 0 ) 

~ Psmlt(T > T 0 ) — Pgmlt (T > T 0 ) + Pjhom(T < T oa ) - Pimlt(T > T a ) + P IDW (T 00 < T < T 0 ) . 

(A4a) 

I-/ Q $ ^ ^ y 

The term ( — --—( — - + — + w — ) ) is omitted in the derivation of (A4) since it is 


r dt dx 


dz 


much smaller than the other terms in (A4). Based on Li et al. (2002c), Pi 8 in a tropical, 
deep convective regime can be simplified as 


14 


"O 11^0 


P„ - P^-PmltoV > < Z - .)- W r > TJ-P^ t (T > TJ. 


(A4fc) 


Taking x (yi2) - = -- V — , the tendency equation for the vertical component 

to p to 

£ dd, . 

of the CPV (P 2 = -=- — L ) can be expressed by 
p to 


ap , ap a p m ? a«afl awao,. lafl.as x / £ a ^0 


at 


dx 

d ( Q R e e) 

a* T ' 


dz 


(A5) 
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Figure Captions 

Fig. 1 Temporal and vertical distribution of vertical velocity (a), zonal wind (b), and time 
series of sea surface temperature (c) observed and derived from TOGA COARE 
for the 10-day period. Upward motion in (a) and westerly wind in (b) are shaded. 
Units of vertical velocity, zonal wind, and sea surface temperature are cm s' 1 , m s' 
l , and °C, respectively. 

Fig. 2 Temporal evolution and zonal distribution of surface rain rates (mm h' 1 ) from 19- 
29 December 1992. 

Fig. 3 Zonal-vertical (X-Z) cross sections of stre amlin es and sum of the mixing ratios of 
cloud hydrometeors (background shading) at (a) 1800 LST, (b) 2100 LST 21, and 
(c) 0000 LST 22 December 1992. 

d0 

Fig. 4 Zonal-vertical (X-Z) cross sections of (a) vorticity (10' 3 s" 1 ), (b) - — - (10' 3 Km' 1 ), 

dz 

df) 

and (c) -—-s- (2x 10^ Km' 1 ) at 2100 LST 21 December 1992. 
dx 

Fig. 5 Time series of zonally-averaged, mass-integrated zonal ([/JJ, lO^Ks' 1 , dot line) 
and vertical ([ ], 10' 3 Ks _I , sobd line) components of the CVY, and sum of 

mixing ratios of cloud hydrometeors ([ g c + q r + q i + q s + q g ], mm, dashed line) 
during the 10-day integration. The plotting scales of [P x ] and [P z ] and 
[q c + q r + q t + q s +q g ] are -0.15-0.15 Ks' 1 , 0-3x 10' 3 Ks _1 , and 0-3 mm, 
respectively. 

Fig. 6 Time series of four contributing processes in the tendency equation for [P z ] in (7) 
during the 10-day integration. Black, green, red, blue, and orange lines denote the 
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tendency of [^], [Term A], [TermB], [Term M], and [TermR], respectively. 
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Fig. 1 Temporal and vertical distribution of vertical velocity (a), zonal wind (b), and time 
series of sea surface temperature (c) observed and derived from TOGA COAJRE for the 
10-day period. Upward motion in (a) and westerly wind in (b) are shaded. Units of 
vertical velocity, zonal wind, and sea surface temperature are cm s' 1 , ms, and °C, 
respectively. 
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Fig- 5 Time series of zonally averaged, mass-integrated zonal ([P ], lO^Ks' 1 , dot line) 
and vertical {[P z ], 1 0 3 Ks ! , solid line) components of the CW, and sum of the mixing 
ratios of cloud hydrometeors ([ q c + q r + q t + q s + q g ], mm, dashed line) during the 10-day 
integration. The plotting scales of [ PJ and [ P z ] and [q c + q r + q i + q s +q g ] are -0. 1 5- 
0.15 Ks" 1 , 0-3 x 10 3 Ks~’, and 0-3 mm, respectively. 
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Popular Summary 


Although dry/moist potential vorticity ( 




) is a useful physical quantity for 


P 

analysis, it cannot be applied to the analysis of 2D simulations. A new potential 


meteorological 
vorticity vector 


— (cloud potential vorticity, CPV) is introduced in this study to analyze 2D cloud-resolving 

P 

simulation data associated with 2D tropical convection. The cloud model is forced by the vertical 
velocity, zonal wind, horizontal advection, and sea surface temperature obtained from the TOGA 
COARE, and is integrated for a selected 10-day period. The CPV has zonal and vertical components 
in the 2D x-z frame. Analysis of zonally-averaged and mass-integrated quantities shows that the 
correlation coefficient between the vertical component of the CPV and the sum of the cloud 
hydrometeor mixing ratios is 0.81, whereas the correlation coefficient between the zonal component 
and the sum of the mixing ratios is only 0.18. This indicates that the vertical component of the CPV 
is closely associated with tropical convection. The tendency equation for the vertical component of 
the CPV is derived and the zonally-averaged and mass-integrated tendency budgets are analyzed. The 
tendency of the vertical component of the CPV is determined by the interaction between the vorticity 
and the zonal gradient of cloud heating. The results demonstrate that the vertical component of the 
CPV is a cloud-linked parameter and can be used to study tropical convection. 


